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ABSTRACT 

We study the dynamical structure of a cooling dominated rotating accretion flow 
around a spinning black hole. We show that non-linear phenomena such as shock waves 
can be studied in terms of only three flow parameters, namely, the specific energy {£ ), 
the specific angular momentum (A) and the accretion rate (m) of the flow. We present 
all possible accretion solutions. We find that a significant region of the parameter 
space in the £ — A plane allows global accretion shock solutions. The effective area of 
the parameter space for which the Rankine-Hugoniot shocks are possible is maximum 
when the flow is dissipation free. It decreases with the increase of cooling effects and 
finally disappears when the cooling is high enough. We show that shock forms further 
away when the black hole is rotating compared to the solution around a Schwarzschild 
black hole with identical flow parameters at a large distance. However, in a normalized 
sense, the flow parameters for which the shocks form around the rotating black holes 
are produced shocks closer to the hole. The location of the shock is also dictated by 
the cooling efficiency in that higher the accretion rate (m), the closer is the shock 
location. We believe that some of the high frequency quasi-periodic oscillations may 
be due to the flows with higher accretion rate around the rotating black holes. 

Key words: accretion, accretion disk - black hole physics-shock waves. 



1 INTRODUCTION 

The accretion onto black holes is thought to be an essen- 
tial process in black hole astrophysics since the gravita- 
tional potential energy released by the infalling matter is 
directly relat ed to the emitted power from the vicinity of 
a black hole IjFrank et al.ll2002T ). In the last three decades, 
a number of works have reported extensive studies of ac- 
cretion disk models around gravitating objects by the sev- 
eral groups both theoretically as well a s through numerical 
simulations (IShakura fc Sunvaevl Il973l ; [Lynden-Belll 1 19781; 
Eggu m et alJ Il9 88: Chakrabarti 1989; No buta fc Hanawal 
1994 1 Yang fcK afatos 1995; Chakrabarti 19 96al;lLu fc. Yuan! 



1997 



Fukumura fc Tsurutall2004l; iNovikov fc Thorndll973l : 



Page fc Thorndll974l ; lHawlev et al.lll984J bl: lGammieHl999h 



In 1980s it was realized that an accreting matter becomes 
transonic around a galactic and extragalactic black hole 
(BH) in order to satisfy the inner boundary conditions at 
the event horizon. A fully general relativistic theoretical 
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work of lChakrabartll l|l996bl ) showed that multi-dimensional 
sub-Keplerian flows exhibit discontinuous shock transitions 
for a significant region of the flow parameter space spanned 
by the specific energy and the specific angular momentum 
in a non-dissipative flow. Recently, it has been pointed out 
that the shock waves formed around a BH could be a very 
useful tool in explaining the observational results. The hot, 
dense post-shock flow could be the natural site for gener- 
ating high en ergy photons jc hakraba rti fc Mandall I2OO6). 
In particular, IChakrabarti fc Mandall (|2006h reported that 
non-thermal electrons may be processed around the shock 
fronts by particle acceleration mechanism that could pro- 
duce emitted radiation up to few MeVs by inverse Comp- 
tonization of non-thermal synchrotron radiations. In addi- 
tion, the observed quasi-periodic oscillations (QPOs) in wide 
frequency range may be originated due to some type of reso- 
nance oscillation of the post-shock flow l|Molteni et al.|[l996l : 
IChakrabarti et al.|[20oi ; lOkuda et all 120071 ), or due to dy- 
nami cal oscillations of the sub-Keplerian flows l|Rvu et alJ 
Il997h . This oscillation is governed by the infall time scale of 
matter in the post-shock region. Since a shock may form very 
close to a spinning BH, it could possibly explain the origin 
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of high frequency QPOs also. Thus, it would be worthwhile 
to study the behavior of the shock waves around a rotating 
BH for more realistic cases by addition of cooling processes 
inside the disk. This is all more important since the post- 
shock pressure which holds the shocks in place, depends on 
the cooling processes. 

So far, major progresses on the theoretical front have 
taken place only when the accretion flow around the rotat- 
ing BH is thin and when the ene rgy dissipation processes 
are inefficient l|Chakrabartilll996bl ). In reality, the accreting 
matter should be affected by the various cooling processes, 
such as bremsstrahlung and synchrotron cooling, most sig- 
nificantly the latter as the magnetic field is expected to be 
present in the disk. Some theoretica l work along this line has 
recently been completed (|Dasll2007r ). However, the effects of 
synchrotron cooling on the dynamical structure of the ac- 
creting flow around a rotating black hole has not yet been 
studied well. Indeed, the study of accretion flows with var- 
ious energy dissipation processes around a compact object 
in the general relativis t ic lim it is very complex. Recently, 
IChakrabarti fc Mondall ()2006h proposed a pseudo-Kerr po- 
tential which describes the space-time geometry around a 
rotating BH quite satisfactorily for the black hole spin pa- 
rameter, au, —1 ^ a-k ^ 0.8. This pseudo-potential, which 
is an extension of the conce pt laid down while writ ing 
the Paczynski-Wiita potential (|Paczvnski fc Wiitalll98oT ) in 
Schwarzschild geometry allows us to study very complex ba- 
haviour of dissipative accreting matter around a rotating 
black hole in a very simple and elegant way. 

In an optically thick and geometrically thin accretion 
flow, cooling processes reduce the bulk energy of the flow 
during accretion while viscosity not only heats the flow 
but also redistributes the angular momentum in the disk. 
The cooling is governed by the density and temperature 
of the flow while viscosity is assumed mainly due to ion 
collisions or its coupling with the mag netic field. A realis- 
tic dis k can have some visco sity, but as IChakrabarti fc Das! 
l|2004 ): IChakrabarti (|l996aT ) already pointed out, standing 
shock exits even in presence of significant viscosity. Thus, in 
this paper we shall explore the properties of shock waves in 
presence of cooling alone. In other words, our focus here is 
to deal with the accretion flows having low angular momen- 
tum upto the marginally bound value. However, we argue 
that basic conclusions regarding our idealized inviscid solu- 
tion would be qualitatively similar to the viscous solution. 
We wish to report viscous solutions elsewhere in the future. 

The suitability of the pseudo-Kerr potential that 
we employ has been well tested for both particle dy- 
namics and fluid dynamics ( Chakrab arti fc Mondall l200rj ; 
iMondal fc Chakrabartil 120061 ') which allows us to use this 
for our present study with sufficient confidence. We iden- 
tify all possible accretion solution topologies including the 
ones which may allow standing shock transitions. In this 
work, we only focus on thin, non-dissipative shock waves. 
We separate various regions in the parameter space accord- 
ing to the solution topologies. We identify a special region 
where standing sho ck conditions (Rankine-H ugoniot shock 
conditions [RHCs]. lLandau fc LifshitzHl959r ) are not satis- 
fied but accretion solution continues to have three critical 
points. Sol utions of t h is kin d generally exhibit oscillating 
shock as in lRvu et al.l (|l997T ) when Schwarzschild geometry 
was used. In the present paper, we ignore bremsstrahlung 



cooling process si nce it is very inefficient compared to the 
synchrotron loss l|Chattopadhvav fc Chakrab arti 200fJ), es- 
pecially in flows around stellar mass black holes. We also do 
not consider Compton cooling effect. These will be carried 
out in future. 

In the next Section, the governing equations are pre- 
sented. In §3 and §4, the accretion solutions are shown. We 
discuss the shock properties in §5. In the final Section, we 
make concluding remarks. 



2 THE GOVERNING EQUATIONS 

In this Section, we present a set of relevant equations that 
governs the dynamical structure of a stationary, thin, vis- 
cous, axisymmetric accretion flow onto a rotating black hole. 
In addition, we assume an ideal equation of state for accret- 
ing gas. The equations are written on the equatorial plane of 
the accretion disk in the frame of the so -called zero angular 
momentu m observers, and are given by jchakrabarti 1996a; 
iDasl 120071 ), 

(i) the radial momentum equation : 



^du + ldP d$ e 
dx p dx dx 



0. 



(ii) the mass flux conservation equation : 

M = T,ux, 

(iii) the angular momentum conservation equation 
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(2) 



(3) 



(4) 



dx Ex dx 
and (iv) the entropy generation equation 

ZuT^=Q+ -Q-. 
ax 

The variables u, p, P denote the radial velocity, density and 
isotropic pressure of the flow respectively. In the present 
work, the effect of gravity is approximated by the pseudo - 
Ker r potential introduced by IChakrabarti fc Mondall (|2006l ) 
and lMondal fc Chakrabarti (2006). The expression of effec- 
tive pseudo-Kerr potential is given by, 



+ - 4$i$ 3 
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Here, x and R represents the 



where, $i 
1 3_ 

R—xq ' x ' 2x' 2 

cylindrical and spherical radial distance from the black hole 
when the black hole itself is considered to be located at 
the origin of the cylindrical coordinate system. The spe- 
cific angular momentum of the flow (which is equivalent 
to —Ufp/ut in a general relativistic flow, tt M being the flow 
velocity components) is denoted by A. Furthermore, xo = 
(0.04 + 0.97a ft + 0.085a! )/2, lu = 2a k /{x 3 + a\x + 2a\) and 
a 2 — (x 2 — 2x + a%)j(x 2 + a\ + 2a\jx). Here, a k represents 
the BH rotation parameter defined as the angular momen- 
tum of the black hole per unit mass. The subscript 'e' stands 
for the quantities calculated at the disk equatorial plan. The 
mass accretion rate is denoted by M and £ and W x </> are the 
vertical averaged density and the viscous stress of the flow. 
The terms s, T, Q + and Q~ represent the specific entropy, 
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the local temperature, the energy gained and the energy lost 
by the flow, respectively. In our model, the accreting flow is 
assumed to be in hydrostatic equilibrium in the vertical di- 
rection of the flow motion. The half thickness h(x) of the 
disk is then obtained by equating the vertical component of 
the gravitational force and the pressure gradient force: 



h(x) 



(6) 



where, <E> R = (■1^1 



and 7 is the adiabatic index of 



' 34^ 

k SRj z<<x 

the flow. Here, z is the vertical component of the distance 
in the cylindrical co-ordinate system and R = y/x 2 + z 2 . 
The a diabati c sound speed is denoted by a and defined as 
a = V "fP/ P- I n t ne weak viscosity limit, W x $ is negligible 
and the angular momentum distribution remains uniform 
throughout the disk. 

The flow equations have been made dimensionless by 
considering unit of length, time and the mass as GMbh /c 2 , 
GMbh/c? and Mbh respectively where, G is the gravita- 
tional constant, Mbh is the mass of the black hole and c 
is the velocity of light respectively. Henceforth, all the flow 
variables are expressed in geometrical units. 

We simplify the entropy equation (Eq. 4) for an inviscid 
flow (Q+ -> 0) as: 



7' 



1 dP_ 

p dx 



7P dp 
p 2 dx 



A. 



(4a) 



where, the term A(= Q~ / ph) represents the energy lost per 
unit gram by the flow. In this work, we only focus on syn- 
chrotron cooling. The strength and the morphology of the 
disk magnetic field are still largely unknown. This motivated 
us to consider only the stochastic magnetic field inside the 
accretion disk and this magnetic field is thought to be at the 
most in equipartition with the accretion plasma. We define 
a control parameter, /3, as the ratio of the thermal pres- 
sure and the magnetic pressure of the flow which acts as a 
measure of equipartition and is given by, 



8TT P k B T p 

B 2 pm p ' 



(7) 



where, B denotes the magnetic field strength, ks is the 
Boltzmann constant, p is the mean molecular weight 
and m p is the mass of the proton respectively. The 
equipartition parameter /3 > 1 in general, ensures 
that magn etic fields remain confined with the accret- 
ing plasma (|Mandal fc Chakra barti 2005]). The synchrotron 
emissivity for the stochas t ic magne tic field is given by 
l|Shapiro fc Teukolskvlll983l ; lDasll2007l ). 



A = 



Sa 5 




with 



S = 15.36 x 10 1 



2 4 
mp e 



/3ml 7 5 / 2 GM Q c 



(8) 



(9) 



where, e and m e represent charge and mass of an electron 
respectively. Here, m is the accretion rate in units of Edding- 
ton rate that regulates the effi ciency of cooling. Following 
iMondal fc Chakrabartil (|2006l ). we use modified polytropic 
index [n = (7 — 1) _1 ] relation as n — > n + (0.3 — 0.2afc) and 
choose (3 = 10 throughout the paper, until otherwise stated. 



3 SONIC POINT ANALYSIS 

We solve Eqs.(l-3, 4 a) following the s tandard method of 
sonic point analysis (jChakrabartl 1989). We calculate the 
radial velocity gradient as: 

du _ N 
dx~~D' 

where, the numerator N is given by, 



(10) 
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2 (7 + 1) d$ e 
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c 3 2(7 — 1) (7 — 1) dx (7 — 1)$^. dx 

(10a) 

and the denominator D is given by, 
2a 2 u 



D 



(7 + 1) 3 
-u . 



(7-1) (7-1) 
The gradient of sound speed is obtained as: 



da /a 71A du 3a 7 d$ e 
dx \u a J dx 2x a dx 



2& R dx 



(106) 



(11) 



Since the matter is accreting onto the BH smoothly except at 
the shock location, the radial velocity gradient must always 
be real and finite. However, eq. (10b) shows that there may 
be some points where denominator (D) vanishes. This indi- 
cates that the numerator (N) must also vanish there to keep 
du/dr finite. These special points where both the numerator 
(N) and denominator (D) vanish simultaneously are called 
critical points or sonic points. Setting D = 0, one can easily 
obtain the expression for the Mach number (M = u/a) at 
the sonic point as, 



M(x c 



7 + 1 



(12) 



We obtain an algebraic equation for the sound speed 
(a c ) by using another sonic point condition N = which is 
given by, 



F(£ c , A c , m) = Aal +Ba c +C = 0, 



(13) 



where, 



A- 



B = 



5(7-1) 



1 3 

x 



& R dx 



and 



C 



(7 + l)M 



2 jge 

dx 



The subscript 'c' denotes the quantities computed at the 
sonic point. We calculate sou nd speed at the sonic point b y 
solving Eq. (13) analytically l|Abramowitz fc Stegunlll970l ). 
In general, a dissipative accretion flow may have multiple 
sonic points depending on the flow parameters. The nature 
of the sonic point is dictated by the sign and the numerical 
value of the velocity gradients at the sonic point. In reality 
du/dr may possesse two real values at the sonic point: one 
is for accretion branch and the other corresponds to the 
wind branch. When both the values of du/dr are real and 
of opposite signs, the sonic point is referred to as saddle 
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Figure 1. Variation of Mach number (M = u/a) with logarithmic 
radial distance. Flow parameters are Xi n = 3.3753, A = 3.04, 
m = 0.0025 and aj. = 0.5. Standing shock forms at x s = 48.57. 

type. A transonic flow generally passes through the saddle 
type sonic points only and for a shock, the flow crosses two 
saddle type sonic points, one before the shock and the other, 
after. The closest one from the BH horizon is called the inner 
sonic point and the furthest one is known as the outer sonic 
point. If the derivatives are real and of same sign, the sonic 
point is nodal type. When both the derivatives are complex, 
the sonic point is of spiral type which is unphysical as no 
physical solution can pass through it. 

4 ACCRETION SOLUTION 

In order to obtain a complete accretion solution, we choose 
the inner sonic point location (xi n ) and th e angular m omen- 
tum (A) of the flow as input parameters l |Dasll2007h . From 
the inner sonic point, we integrate inward up to the BH hori- 
zon and outward till the outer edge and combine them to 
obtain a global transonic solution. 

4.1 Shock Solution 

In Fig 1, we present a global solution with a standing shock 
around a rotating black hole. The variation of the Mach 
number with the logarithmic radial distance is plotted. The 
flow parameters are Xi n — 3.3753, A = 3.04, m = 0.0025 and 
at = 0.5 respectively. The figure consists of two solutions. 
The one passing through the outer sonic point (O) connects 
the black hole horizon and the outer edge of the disk. The 
other passing through inner sonic point (I) is closed and is 
connected with the BH horizon only. Arrows indicate the di- 
rection of the flow towards the BH. Matter starts accreting 
from the outer edge of the disk with a negligible velocity. 




X 

Figure 2. Parameter space for multiple sonic points. The effec- 
tive region bounded by the solid curve are for cooling free accre- 
tion flow. Regions under dotted and dashed curves are obtained 
for m = 0.0025 and 0.0125 respectively. Effective area of param- 
eter space reduces with the increase of accretion rate (m) as it 
enhances the cooling efficiency. 

As the flow accretes towards the BH, the flow gains its ra- 
dial velocity due to the attraction of the strong gravity. The 
flow becomes super-sonic after crossing the outer sonic point 
(O) and continues to acc rete towards the BH. The RHCs 
IlLandau fc Lifshitd[l959l) in turn allows the flow variables 
to make a discontinuous jump in the sub-sonic branch. This 
is indicated by a dotted vertical line and commonly known as 
the shock transition. In the post-shock region, the flow mo- 
mentarily slows down and subsequently picks up the radial 
velocity and enters into the BH supersonically after crossing 
the inner sonic point (I). In this particular case, the shock 
conditions are satisfied at x s — 48.57. The entropy generated 
at the shock is eventually advected towards the BH to allow 
the flow to pass through the inner sonic points. In this work, 
the shocks are considered to be thin and non-dissipative. 



4.2 Parameter space for multiple sonic points 

So far, we have seen that the shock wave connects two so- 
lution branches — one passing through the outer sonic point 
and the other passing through the inner sonic point. In par- 
ticular, the flow with multiple saddle type sonic points may 
undergo shock transitions. More importantly, for shock for- 
mation in dissipative accretion flows, the solution pa ssing 
through the inner sonic point has to be spiraling in (IDasI 
l2007h . Therefore, it would be useful to study the parameter 
space for accretion flows having spiral in solution passing 
through the inner sonic point. In Fig. 2, we show the classi- 
fications of parameter space as a function of accretion rate 
(rh) in the £i„ — A plane. Here, £i„ denotes the energy of the 
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Figure 3. Mach number variation with logarithmic radial dis- 
tances. Flows are injected sub-sonically from the outer edge 
%inj = 300 with identical energy Ei n j = 1.003163 and angular 
momentum A = 3.0. Different accretion rates (rh) are used. Solid 
curve represents a solution including the shock wave (x a = 26.03) 
for cooling free accretion flow. Other solutions are for [rh, x s ]= 
[0.0125, 19.38] (dotted), [0.025, 15.85] (dashed) and [0.0375, 13.45] 
(dot-dashed). As rh is increased, shock front precedes towards BH. 

flow at the inner sonic point (£,„). The BH rotation param- 
eter is chosen as at = 0.5. The region bounded by the solid 
curve is obtained for non-dissipative accretion flows. The 
regions under the dotted and dashed curve is obtained for 
higher accretion rates rh = 0.0025 and 0.0125 respectively. 
As accretion rate (rh) is increased, the parameter space for 
multiple sonic points shrinks. This indicates that nature of 
sonic points changes (from saddle to spiral) for flow with 
identical input parameters for increasing rh. 



5 SHOCK PROPERTIES 
5.1 Shock Dynamics 

In Fig. 3, we present the variation of shock locations with 
the accretion rate (rh). Logarithmic radial distance is varied 
along the horizontal axis and Mach number is plotted along 
vertical axis. The vertical lines represent the shock locations. 
Matter with identical outer boundary conditions is injected 
sub-sonically from the outer edge of the disk Xi„j = 300 on 
to a rotating BH with rotation parameter = 0.5. The 
local energy of the flow at x in j is £ in j = £(x in j) = 1.003163 
(including the rest mass) and the angular momentum is 
A = 3.0. The sub-sonic flow crosses the outer sonic point 
to become super-sonic and makes a shock transition to the 
sub-sonic branch. The solid vertical lines represent the shock 
location (x s — 26.03) for non-dissipative flow. As cooling 
is incorporated, shock front moves forward. In the post- 
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Figure 4. (a) Variation of the shock location with accretion rate 
(fa). Flows arc injected with the same energy and angular mo- 
mentum from the outer edge. Small dashed, big dashed, small-big 
dashed, dot-dashed, dotted and solid curves are drawn for angu- 
lar momentum A = 2.94, 2.96, 2.98, 3.0, 3.02 and 3.04 respectively, 
(b) Variation of compression ratio (R = S+/S_) with accretion 
rate for the same set of parameters as in (a), (c) Variation of 
shock strength (0 = iW_/M_|_) with the accretion rate for the 
same set of parameters as in (a). Subscripts "+" and "-" denote 
quantities before and after the shock. 

shock region, cooling is more effective compared to the pre- 
shock flow as the density as well as the temperature are 
very high in this region due to compression. Cooling reduces 
the post-shock pressure causing the shock front to move 
inward to maintain pressure balance across it. For higher 
rh, shock front moves further inward. The dotted, dashed 
and dot-dashed vertical lines represent the shock locations 
x s = 19.38, 15.85 and 13.45 for rh = 0.0125, 0.025 and 0.0375 
respectively. 

In Fig. 4a, we show the variation of the shock location 
as a function of accretion rate (rh) for a set of angular mo- 
mentum (A). In this particular figure, the flow is injected 
from the outer edge of the disk (xi n j = 300). The BH ro- 
tation parameter is chosen as = 0.5. The angular mo- 
mentum of the flow is varied from A = 2.94 (small dashed) 
to 3.04 (solid) with an interval AA = 0.02. At the injection 
point, the corresponding local energies of the flow (from bot- 
tom to top) are £ inj = 1.003174 (small dashed), 1.003170 
(big dashed), 1.003167 (small-big dashed), 1.003163 (dot- 
dashed), 1.003160 (dotted), and 1.003156 (solid) respec- 
tively. For a given accretion rate (rh), the shock forms fur- 
ther out for flows with higher angular momentum. Here, the 
larger angular momentum increases the centrifugal pressure 
which pushes the shock front outside. Conversely, for a given 
angular momentum, the shock location decreases with the 
increase of the accretion rate as cooling reduces the post- 
shock thermal pressure. Figure 4a shows that the standing 
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shocks are formed for a wide range of accretion rate (rh). 
In each angular momentum, the standing shocks disappear 
beyond a critical value of accretion rate (rh c ) as the RHCs 
are not satisfied here. Non-steady shocks may still exit, but 
an investigation of such phenomena is beyond the scope of 
the present work. As A increases, the critical accretion rate 
first increases, becomes maximum at some A (= 2.96, in this 
particular case) and then decreases. This clearly indicates 
that the parameter space for the standing shock shrinks in 
both the lower and higher angular momentum sides with the 
increase of the accretion rate. 

One of the important components in accretion disk 
physics is to study the density profile of matter since the 
cooling efficiency as well as the emitted radiation directly 
depends on it. We compute the compression ratio R defined 
as the ratio of vertically averaged post-shock to pre-shock 
density and plot it in Fig. 4b as a function of accretion rate 
(rh) for the same set of flow parameters as in Fig. 4a. For 
a given angular momentum (A), the compression ratio in- 
creases monotonically with higher rh. As rh increases, post- 
shock flow becomes more compressed to provide required 
pressure for holding the shock. In addition, for a given rh, 
higher angular momentum flow feels less compression in the 
post-shock region as centrifugal pressure resists the flow to 
accrete. Note that, for each angular momentum, there is a 
cut-off at a critical accretion rate limit as standing shock 
conditions are not satisfied there. 

It is useful to study the another shock property called 
shock strength O (defined as the ratio of pre-shock to post- 
shock Mach number of the flow) as it is directly related to 
the temperature jump at the shock. In Fig. 4c, we show the 
variation of shock strength as a function of accretion rate 
(rh) for flows with identical input parameters as in Fig. 4a. 
For a given angular momentum (A), the strength of the shock 
is the weakest in the dissipation-free limit and it becomes 
stronger as accretion rate (rh) is increased. Thus, a higher 
cooling causes the post-shock flow to be hotter and radia- 
tions emitted from this regi on are expected to be harder. A 
similar result is reported bv lMandal fc Chakrabartil (2005). 
This clearly indicates that the observed spectra of the BH 
would strongly depend on the cooling. 

An important part of understanding a cooling domi- 
nated accreting flow around a rotating BH is to study the 
shock properties as a function of BH rotation parameter ak- 
in Fig. 5a, we plot the variation of shock location as a func- 
tion of afe. In this particular Figure, we inject matter from 
the outer edge x% n j = 200 and the accretion rate is consid- 
ered to be rh = 0.0025. Solid, dotted and dashed curves are 
obtained for flows with angular momentum A = 2.96, 3.05 
and 3.14 respectively. The corresponding energies at the in- 
jection point are £ inj = 1.003456, 1.003442 and 1.003422 
respectively. Notice that, for a given A, shocks form for a 
particular range of ak- As at increased, shock recedes from 
the BH horizon. Moreover, shocks exist around a weakly ro- 
tating BH when the flow angular momentum is relatively 
higher and vice versa. This phenomenon is directly related 
to the spin-orbit coupling term in the Kerr geometry. In fact, 
since both the marginally bound and the marginally stable 
angular momenta (as well as their difference) go down when 
the Kerr parameter is increased, the relevant parameter re- 
gion when the shocks form also goes down as ak is increased. 
In general, however, the shock location is generally small for 
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Figure 5. Variation of (a) shock location, (b) compression ratio 
and (c) shock strength with BH rotation parameter a^. See text 
for more details. 

higher ak, as statistically the flow with a smaller angular 
momentum (and therefore, the lesser centrifugal force) is 
accreted in a rapidly spinning black hole. Thus, for instance 
if we compared the shock locations having A = A m3 for all 
cases, the shock location for a rotating black hole would be 
closer for higher ak- 
in Fig. 5b, we show the variation of the compression 
ratio as a function of BH rotation parameter at for the flows 
with input parameters same as Fig 5a. The solid, dotted 
and dashed curves are obtained for A = 2.96, 3.05 and 3.14 
respectively. The compression ratio R decreases with the 
increase of ak for flow with identical A. In Fig. 5c, we plot 
the variation of shock strength O with ak for flow with input 
parameters as in Fig. 5a. We obtain a similar variation of G 
with a^ as in Fig. 5b. 

5.2 Parameter Space for Shock Formation 

In Fig. 6, we identify the region of the parameter space that 
allows the formation of the standing shocks. The BH ro- 
tation parameter is considered to be ak = 0.5. The region 
bounded by the solid curve is obtained for non-dissipative 
accretion flow. As accretion rate is enhanced, the effective 
region of parameter space for standing shocks shrinks in 
both the lower and higher angular momentum sides. Due 
to the cooling effect, the flow loses its energy as it accretes 
and therefore, the parameter space is shifted to the lower 
energy domain for higher cooling. The regions under dot- 
ted, dashed and dot-dashed curves are obtained for accretion 
rates rh = 0.0025, 0.005 and 0.01 respectively. It is clear that 
the standing shocks do not exist beyond a critical accretion 
rate when the synchrotron cooling is present. 

In Fig. 7, we classified the entire parameter space 
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Figure 6. Variation of the effective region of parameter space 
which forms standing shocks as a function of the accretion rate 
(rii). See the text for more details. 
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Figure 7. Classification of parameter space according to the var- 
ious solution topologies of BH accretion solution. See text for 
details. 

spanned by (£i„,X) according to the nature of solution 
topologies. As an example, we consider = 0.5 and m = 
0.0025. We separate the parameter space into six regions 
marked by S, OS, O, I, CI and N. The dot-dashed line rep- 
resents the rest mass energy of the flow. At the bottom left 
of the parameter space, we plot solution topologies in the 
small boxes. In each box, Mach number of the flow is plot- 
ted against the logarithmic radial distance. Each of these 



solutions are marked and drawn using the parameters from 
the corresponding region of the parameter space. The direc- 
tion of the accreting flow is indicated by the arrow. (The 
solutions without arrows are relevant for winds, discussions 
on which are beyond the scope of this paper.) The solutions 
from the regions marked 'S' and 'OS' have two X type sonic 
points and the entropy at the inner sonic point is higher 
than that at the outer sonic point. Flows from 'S' suffer a 
standing shock transition as RHCs are satisfied. However, 
a solution from the region 'OS' does not pass through the 
standing shock as RHCs are not satisfied here. A flow with 
parameters from this region is unstable and causes periodic 
variation of emergent radiation from the inner part of the 
disk as it tries to make a shock transition but fails to do 
so. This is kno wn from th e num erical simulations of non- 
dissipative flows iRvu et al.l (|1997T ) and we anticipate that a 
similar behaviour would be seen in this case as well. The so- 
lutions from the region T possess only the inner sonic point 
and the accreting solutions straight away pass through it 
before entering into the BH. A solution from the region 'O' 
has only one outer sonic point. The solution from region 'CI' 
has two sonic points — one 'X' type and other 'spiral' type. 
Solutions of this kind does not extend to the outer edge of 
the disk to produce a complete global solution a nd therefore, 
become s unstable. It has been pointed out by IChakrabartil 
(|l996aT ) that inclusion of viscosity should open up the topol- 
ogy to allow the flow to reach a larger distance to join with 
a sub-sonic Keplerian flow. The region marked 'N' is the 
forbidden region for a transonic flow solution. 



6 CONCLUDING REMARKS 

In this paper, we have studied the properties of cooling dom- 
inated accretion flow around a rotating black hole by solv- 
ing a set of equations that regulate the dynamical structure 
of the flow. A special consideration is given to synchrotron 
cooling that strongly affects the disk properties as well as 
the emitted spectrum and luminosity that are observed. 

We obtain the global accretion solutions with and with- 
out the shocks in terms of a few flow parameters, namely, 
energy, angular momentum, black hole rotation parameter 
and the accretion rate, which effectively acts as the cooling 
parameter. We find that the accreting matter experiences a 
centrifugal force which acts as a barrier, inducing a shock 
formation. We show that the global shocked accretion solu- 
tion can be obtained for a significant region of the parameter 
space even when the cooling is significant. Using a conven- 
tional accretion disk model we expect the accretion to take 
place when the angular momentum is close to the marginally 
stable value. Our calculation shows that the region is actu- 
ally broader, in terms of both the angular momentum and 
energy. 

The discussion regarding the natur e of the sonic point 
has been reported in many oc casions l|Chakrabarti fc Das! 
|2004| ; iDas fc Chakrabarti|[2004r i. However, a detailed anal- 
ysis was not presented before for a cooling dominated flow 
around a rotating black hole. Our present work suggests that 
a large region of the parameter space provides a stable sad- 
dle type sonic point. In Fig. 2, we demonstrated that the 
parameter space for the stable saddle type sonic point is 
gradually reduced with the increase of cooling efficiency. 
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We show that the standing shocks form closer to a spin- 
ning black hole as the accretion rate is enhanced. At the post 
shock region, the density and the temperature is relatively 
high compared to the pre-shock flow and thus cooling is more 
efficient there. For a higher cooling, the post-shock matter 
cools faster reducing the thermal pressure drastically. This 
forces the shocks to move inward to maintain pressure bal- 
ance across them. 

One of the aims of the present work was to study the 
effect of black hole rotation parameter on the dynamical 
structure of cooling dominated global solutions. We find 
that for flows with identical outer boundary condition (e. 
g., same energy and angular momentum at the outer edge) 
shock recedes from the black hole horizon with the increase 
of black hole rotation parameter (ctfc). However if we choose 
the relevant angular momentum for each case, such as the 
marginally stable angular momentum the shock location 
moves in with the increase of au- The range of at for which 
the stationary shocks are formed is restricted for a flow 
of given angular momentum. Shocks are possible around a 
rapidly rotating black hole when the flow angular momen- 
tum is relatively low. Since that produces a very low cen- 
trifugal pressure, the shock can form very close to the black 
hole for a rapidly spinning black hole. 

We identify the region of parameter space for the forma- 
tion of a standing shock. We find that the effective region of 
the parameter space for the stationary shock shrinks when 
the accretion rate is enhanced. This suggests that the possi- 
bility of shock formation decreases for higher accretion rate. 
In addition, we also separate a region where the Rankine- 
Hugoniot relation is not satisfied. In the context of invidcid 
flows, it has been observed that the flow parameters from 
such a region give rise to oscillating shocks l|Rvu et al.ll 19971 ) . 
The reason is that the higher entropy at the inner sonic 
point forces the flow to pass through it by generating extra 
entropy at the shock. But since RHCs are not satisfied the 
shock can not settle itself at a given location. Thus the cause 
of oscillation is sufficiently generic and we suspect that ex- 
actly the same thing will happen in the present case. Most 
importantly, since rotating black holes may have shocks very 
close to the horizon, the frequencies of such oscillations are 
expected to be higher. 

Our present findings suggest that shocks, standing or os- 
cillating, do form around the spinning black holes and it may 
be an essential ingredient sinc e shocks could successfully ex - 
plain the observed stationary (Chakrab arti fc Mandalll2006l ) 
as well as time dependent beh aviour of the radiations 
from the black h ole candidates (jChakrabarti et al . 2004; 
lOkuda et af]|2007h . We demonstrated that the shocks form 
closer to the black hole as cooling is increased. This will en- 
hance the QPO frequency as it is proportional to the infall 
time scale (|Chakrabarti fc Manick am 2000; iMolteni et al.l 
1996) and thus the QPO frequency may vary in a wide 
range starting from mHz to KHz depending on the accre- 
tion rate. This understanding also generally agrees with the 
observational results. Recent reporting of the outbursts of 
GRO 1655-40 showed a clear evidence of the QPO frequen- 
cies increasing monotoni cally from about 90mHz to 17Hz 
jChakrabarti et al.l [iooBh in a matter of 15 days and this 
could be easily fitted using the shock propagating at a con- 
stant velocity (Chakrabarti et al. 2005). 

The formalism presented here does not include out- 



flows/jets which may be generated from the inner part of 
the disk as a result of deflection of inflowin g matter due to 
excess thermal pressure at the shock front (IDas et alj|200ll ; 
IChattopadhvav fc DaslboOTl ; IDas fc Chattopadhvavl l2008h . 
Since the outflows/jets are ejected evacuating the inner part 
of the disk, it will necessarily reduce the post-shock pressure 
and therefore, shock front has to move in to retain the pres- 
sure balance. This suggests that the result should be affected 
if the accretion-ejection mechanism is considered together. 
We plan to consider this study in a future work and it will 
be reported elsewhere. 

We have approximated the effect of general relativity 
using the pseudo-Kerr gravitational potential. This pseudo- 
Kerr potential has been successfully tested to retain most, if 
not all, of the salient features of the flows in a Kerr metric. 
The use of this approach allows us to find out the non-linear 
shock solutions in a curved space-time geometry in a sim- 
pler way. We believe that our basic results would be quali- 
tatively the same with fully general relativistic calculations, 
especially for au < 0.8 for which the pseudo-Kerr potential 
was found to be satisfactory. 
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